Calorimetric Evidence for a Fulde— Ferrell— Larkin— Ovchinnikov Superconducting State 
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The specific heat of the layered organic superconductor k-(BEDT-TTF) 2 Cu(NCS) 2 , where BEDT- 
TTF is bisethylenedithio-tetrathiafulvalene, has been studied in magnetic fields up to 28 T applied 
perpendicular and parallel to the superconducting layers. In parallel fields above 21 T, the su- 
perconducting transition becomes first order, which signals that the Pauli-limiting field is reached. 
Instead of saturating at this field value, the upper critical field increases sharply and a second first- 
order transition line appears within the superconducting phase. Our results give strong evidence 
that the phase, which separates the homogeneous superconducting state from the normal state is a 
realization of a Fulde-Ferrell-Larkin-Ovchinnikov state. 

PACS numbers: 74.70.Kn, 75.40.Cx, 74.25.Bt, 65.40.Ba 



The influence of high magnetic fields on the supercon- 
ducting state is of high technological and fundamental 
relevance. Applied aspects such as high superconduct- 
ing current densities as well as basic questions on the 
nature of the superconducting state are of special inter- 
est. In type-II spin-singlet superconductors, a magnetic 
field destroys superconductivity in two distinct ways: i.e., 
by orbital and Pauli-paramagnetic pair-breaking effects. 
In general, both effects limit the maximum upper criti- 
cal field, H C 2, above which the normal-conducting state 
is restored. In most typc-II superconductors the orbital 
effect plays the dominant role. If, however, under cer- 
tain circumstances, the orbital pair-breaking field, H or b, 
is clearly larger than the Pauli paramagnetic limit, Hp, 
the superconductor may enter a special, spatially modu- 
lated, superconducting state at high magnetic field and 
low temperatures. As predicted independently by Fulde 
and Ferrell [l[ and Larkin and Ovchinnikov 0] in 1964, in 
this so-called FFLO (or LOFF) state superconductivity 
can survive even above the Pauli limit by "sacrificing" 
parts of the material volume to the normal state. 

In the FFLO state, the Zeeman-split Fermi surfaces 
allow Cooper pairing only with a finite center-of-mass 
momentum q resulting in an oscillating part of the order 
parameter in real space with wavelength of the order of 
the coherence length, £. For the FFLO state to occur the 
so-called Maki parameter a — y/2H or b/Hp, should be 
larger than 1.8 [3] and the superconductor needs to be 
in the clean limit with a mean-free path, £ , much larger 
than £. Not many superconductors fulfil these conditions 
and early searches failed to observe the FFLO state. In 
the 1990s, some reports claiming evidence for the FFLO 
state appeared that later had to be revised. For more 
details see Q and the recent reviews @, 0] • 

Only recently, solid thermodynamic evidence for the 



existence of the FFLO state has been put forward for the 
heavy- fermion compound CeCoIns [1, [9| . Besides heavy- 
fermion superconductors the quasi-two-dimensional (2D) 
organic superconductors have been suggested as good 
candidates for exhibiting the FFLO state @, 0, E11H • 
For these the orbital pair breaking can be greatly sup- 
pressed when applying the magnetic field parallel to the 



highly conducting layers [lfj, lll|, [12| . Indeed some signs 



but no thermodynamic proof, for an FFLO state in 2D 
organic superconductors have been reported HI, 14, HI- 
For A-(BETS) 2 GaCLt a kink in the thermal conductiv- 
ity [3 and for A-(BETS) 2 FeCl 4 dip structures in the 
resistance (THJl suggested the existence of a FFLO state 
(BETS is biscthylcncdithio-tetrasclcnafulvalcnc). For 
k-(BEDT-TTF) 2 Cu(NCS) 2 the existence of the FFLO 
phase had been inferred from measurements sensitive to a 
loss in vortex stiffness (BEDT-TTF is bisethylenedithio- 
tetrathiafulvalene) The observed features are, how- 
ever, rather broad anomalies that do not fit with our 
specific-heat results discussed below and are, therefore, 
most probably not related to the FFLO state. 

In this Letter, we present clear thermodynamic evi- 
dence that for k-(BEDT-TTF) 2 Cu(NCS) 2 a narrow in- 
termediate superconducting state, most probably an 
FFLO state, evolves at high magnetic fields applied par- 
allel to the layers. At this Pauli-limited field region the 
slope of the upper-critical-field line increases sharply and 
a first-order transition appears below i? c2 . 

Single crystals of k-(BEDT-TTF) 2 Cu(NCS) 2 were 
grown by the standard electrochemical-oxidation method 
as described elsewhere [l6j]. The specific heat was mea- 
sured by use of a miniaturized relaxation technique 
17, HH in Geneva up to 14 T and in the Grenoble High 
Magnetic Field Laboratory up to 28 T. In both cases, the 
same calorimeter was used without removing the sam- 
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pie between the experiments. The chip resistance and 
the thermal conductance of the leads have been carefully 
calibrated up to 28 T using a capacitance thermometer. 
Each relaxation provides about 1000 data points over a 
temperature interval of 30-40% above the base tempera- 
ture, which has been varied between 1.3 and 12 K. Data 
can be recorded during heating and cooling, which al- 
lows to resolve hysteresis effects close to first-order tran- 
sitions. The merging of the upward and downward relax- 
ation data provides a highly reliable check of the accuracy 
of this method. 

The specific heat is a purely thermodynamic bulk 
quantity. Effects related to flux pinning or vortex insta- 
bilities can thus be excluded [l9(. We precisely aligned 
the sample for parallel-field orientation by slightly turn- 
ing the cryostat in the bore of the resistive coil while 
maximizing the superconducting transition temperature 
in an intermediate field of 8 T. We found that a good 
indicator for a perfect orientation was the absence of the 
first-order vortex-melting transition in fields between 3 
and 14 T which is observable in slightly tilted fields [20I ]. 

In order to determine the particular band-structure pa- 
rameters of the investigated specific-heat sample we mea- 
sured the de Haas-van Alphen (dHvA) effect by use of 
the cantilever-torque technique in a 3 Hc cryostat. 

The inset of Fig.QJa) shows specific- heat data taken in 
and 14 T applied perpendicular to the superconduct- 
ing layers. The specific-heat anomaly, which represents 
only 5% of the total specific heat, is visible in zero field 
at T c = 9.1 K. The curves merge for all fields between 4 
and 14 T which indicates that the upper critical field is 
reached at about 4 T for this field orientation. The 14- 
T data thus represent the normal state and allow us to 
extract the Sommerfcld constant, 7, and to separate the 
phonon background. We obtain 7 = 26(2) mJmol -1 K -2 
and a Debye temperature of about 200 K. These val- 
ues as well as the overall specific heat are in very good 
agreement with earlier results 21, 22, [3]. Subtraction 
of the phonon contribution, i.e., the 14-Tesla data minus 
7T, results in the electronic specific heat, C e , shown for 
fields applied perpendicular and parallel to the layers in 
Fig. Hlja) and [ljb), respectively. From the size of the 
specific-heat jump AC = 0.58 Jmol _1 K -2 we obtain the 
ratio AC/(^T C ) = 2.5(2), which, in accordance with ear- 
lier reports 0, [23[, is larger than the BCS value (1.43) 
and thus proves strong coupling. 

The specific heats look rather different for the two 
field orientations. For perpendicular fields the supercon- 
ducting transition is strongly broadened already in small 
fields due to fluctuations. This can be explained by a 
field-induced finite-size effect in layered superconductors 
24j . Contrary to standard superconductors in which T c 
is lowered in a field but essentially remains sharp, the 
onset of the transition in k-(BEDT-TTF) 2 Cu(NCS) 2 is 
hardly influenced. Instead, each sign of superconduc- 
tivity continuously fades away while reaching _ff c2 . A 




FIG. 1: (color) Temperature dependence of the electronic 
contribution to the specific heat, C e /T, of k-(BEDT- 
TTF)2Cu(NCS)2 in magnetic fields applied perpendicular (a) 
parallel (b) to the superconducting layers. Data taken in 14 T 
applied perpendicular to the layers (inset) were used to sepa- 
rate the phonon contribution. 



possible scenario may therefore be that in k-(BEDT- 
TTF) 2 Cu(NCS)2 only the phase coherence is lost in mag- 
netic fields, while Cooper pairs still exist above i? c2 . 

For fields applied parallel to the layers we find as well 
a broadening of the superconducting transition, but this 
effect is less pronounced and a clear shift of the specific- 
heat anomaly occurs. This resembles more a standard 
behavior. The effect of fluctuations is nevertheless also 
visible for this field orientation. 

Above ~14 T, T c shifts rather rapidly with increasing 
field and the transition becomes strongly broadened by 
fluctuations. At about 21 T, the anomaly in C sharpens 
and at 21.5 T a spike in C due to the latent heat of a 
first-order transition appears (Fig. [2]). For higher fields 
up to 23 T, two sharp anomalies in C clearly prove the 
existence of an additional thermodynamic phase within 
the superconducting state. The lower transition shows a 
well-resolvable temperature hysteresis (^0.05 K at 22 T), 
while the main superconducting transition develops only 
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FIG. 2: (color) Details of C e /T of k-(BEDT-TTF) 2 Cu(NCS) 2 
in high magnetic fields applied parallel to the superconducting 
layers. 
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FIG. 3: (color) Phase diagram of k-(BEDT-TTF) 2 Cu(NCS) 2 
for fields applied parallel to the superconducting layers. The 
solid line represents the calculated _ff c2 using the known band- 
structure parameters. See the text for details. 



a small hysteresis < 0.02 K (inset of Fig. We were 
able to follow the two transitions up to 23 T, above which 
they left the temperature window of our experiment. 

The extracted magnetic phase diagram with £f c2 and 
the second transition including its hysteresis is shown in 
Fig. [31 At high temperatures, if c2 increases very steeply, 
before it levels off towards saturation at lower T. This 
clearly signals the crossover from an orbital T c reduction 
at low fields towards a Pauli-paramagnetic limitation at 
higher fields. From the initial critical-field slope, H' c2 = 
dH C 2/dT, the orbital- limiting field, H or {,, can be esti- 
mated. Using the T c reduction of about 0.4 K at 8 T, 
i.e., H' c2 = 20 T/K, wc obtain H orb = 0.7H' c2 T c w 130 T, 
much larger than the low-temperature H a i we observe. 

The large value of H rb results from the strongly re- 
duced orbital currents for the field orientation parallel 



to the layers. For that reason the field-induced spin po- 
larization becomes important resulting in the rapid T c 
reduction. For k-(BEDT-TTF) 2 Cu(NCS) 2 , the Pauli- 
limiting field can be determined quite accurately from 
Hp = A /(V2hb) HI, where fiB is the Bohr magneton 
and the superconducting energy gap, Ao, is well known 
from specific-heat studies [IH l23|. Using A /kpT c = 2.4 
[23 | . we obtain Hp = 23 T which agrees very well with 
the observed limitation of H C 2 towards low temperatures 
(Fig. [3|). fndeed, the crossover of the broadened su- 
perconducting transition to a sharp first-order transition 
shows that the Pauli-limiting field is reached at ~21 T. 

Above this field, the H c2 line clearly increases its slope 
and the transition develops a latent heat. Simultane- 
ously, the second transition line appears within the super- 
conducting phase. This strongly suggests the evolution of 
the FFLO state in k-(BEDT-TTF) 2 Cu(NCS) 2 . Indeed, 
this superconductor fulfils all requirements necessary for 
the formation of the FFLO state. First, it is a strongly 
type-II superconductor with a large Ginzburg-Landau 
parameter k of 100 — 200 even for perpendicular fields 
[261 ]. Second, the Maki parameter, a = V2H or b/Hp w 8 
is more than four times larger than required [4| . Finally, 
£ is much smaller than the mean-free path, £. The lat- 
ter was determined for the investigated sample by means 
of dHvA measurements. In accordance with well estab- 
lished literature data [27j], two dHvA frequencies were 
resolved. The lower frequency, F a = 600(f) T, orig- 
inating from a hole orbit, together with the measured 
effective mass, = 3.05(I0)m e , allows to calculate the 
Fermi energy t c p = 23 mcV, where m e is the free-electron 
mass. For the larger breakdown orbit \Fp = 3870(20) T, 
mP c = 6.5(2)TO e ], comprising all electrons at the Fermi 
level, the Fermi energy is e F = 69 meV. The field depen- 
dence of the dHvA amplitude gives the scattering time 
t = 1.9 x I0~ 12 s for the a orbit. The magnetic break- 
down allows the determination of the scattering time for 
the (3 orbit only with much larger error bars. Assuming 
an unchanged t leads to £ = 100 and 115 nm, respec- 
tively. By use of H C 2± = 4 T from our specific-heat data 
[Fig. DJa)], £ f=a 9 nm. Correspondingly, £/£ rts 12 proves 
that our sample is in the clean limit. 

We may now compare our results with the calcu- 
lated phase diagram for magnetic fields parallel to the 
layers. For that, we assume s-wave superconductivity 
[23 . [23l . [28j], although that is not an essential ingredi- 



ent. Detailed studies [U|, [12j showed that the transi- 
tion lines, the orders of the transitions, and the struc- 
tures of the equilibrium states depend very sensitively 
on the electronic structure, i.e., on the Fermi surface, 
the effective masses of the quasiparticles, and their in- 
teractions. We model the system under consideration 
by a stack of 2D superconducting planes with negligi- 
bly small conductivity perpendicular to the planes and 
use the Fermi surface and effective masses derived from 
the dHvA experiments. As the magnetic susccptibil- 
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ity in the normal state is not significantly renormalized 
compared to the expected value for free quasiparticles 
we neglect Landau's spin-dependent interaction parame- 
ter Fq. For a quantitative comparison with experiment 
we have to account for corrections due to the strong 
electron-phonon coupling. Considering the rather small 
ratio (kBT c )/(h(cu)) ~ 0.08 [29[ we anticipate neither 
pronounced anomalies in the variation with T of the en- 
hancement factors nor significant differences in critical- 
field renormalizations corresponding to first and second- 
order transitions, respectively [30l |. As we are mainly 
interested in the low-T behavior we implicitly account 
for strong-coupling corrections by rescaling the weak- 
coupling results with the enhancement factor for the low- 



T energy gap A /(1.76/c B T c ) ~ 1.35 [23J]. The solid 
line in Fig. [3] displays the resulting critical field for the 
second-order transition from the normal to the inhomo- 
geneous FFLO state which agrees well with experiment. 
The anisotropy of the effective masses which differ by a 
factor of ~2 on the different Fermi-surface sheets stabi- 
lizes the inhomogeneous superconducting state and leads 
to the steep upturn in H C 2(T). Following Ref. [3l[ we 
find that the transition is of second order in the vicinity 
of the tricritical point in close analogy to the case of the 
isotropic 2D superconductor with non-interacting quasi- 
particles. For a better description of the low-T H C 2(T) 
data and for estimating the transition from the homoge- 
neous superconducting to the FFLO state more detailed 
information on the electronic parameters is necessary. 

When we compare the phase diagram in Fig. [3] with 
that of CeCoIns [5J] we find clear qualitative differences. 
For CeCoIn 5 , H c2 is less reduced at high fields due to 
Pauli limitation and there is no clear upturn of H C 2 when 
the FFLO state appears. The latter may be due to the 
rather isotropic in-plane effective masses. For the 2D or- 
ganic material, the FFLO phase occupies only a small 
area in the phase diagram with the transition line follow- 
ing closely the H C 2 transition (Fig. [3]) . Another difference 
is found in the nature of the transitions. For the organic 
superconductor both transitions are found to be first or- 
der, whereas for CeCoIn 5 only the H c2 line is first order. 
Both materials have similar values for a and £/£ (see [32[ 
for recent estimates). A possible reason for these dif- 
ferences might be the anisotropy which for CeCoIns is 
much smaller than for k-(BEDT-TTF) 2 Cu(NCS) 2 . For 
the latter the ratio of the Fermi energy with respect to 
the interlayer hopping, ep/t±,is about 3700 (33j. Finally, 
we see no evidence for an FFLO state in the field per- 
pendicular orientation, whereas for CeCoIns the matter 
is still unresolved (HEH- 

In conclusion, we presented high-resolution specific- 
heat data in high magnetic fields that give clear thermo- 
dynamic evidence for the existence of a narrow additional 
superconducting phase in the 2D organic superconductor 
k-(BEDT-TTF) 2 Cu(NCS) 2 for in-plane magnetic fields. 
This phase is most probably a realization of the long- 



time predicted FFLO state. This is supported by the 
paramagnetic limitation of H C 2 , the upturn of i/ c2 when 
the Pauli-limiting field of 21.5 T is reached, and the first- 
order nature of the transitions above this field. The clear 
observation of the FFLO state in a second superconduc- 
tor, besides CeCoIns, with a qualitative different phase 
diagram allows for thorough tests of our fundamental un- 
derstanding of superconductivity at high magnetic fields. 
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